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We summarize and compare recent Molecular Dynamics simulations on the interactions of dipalmi- 
toylphosphatidylcholine (DPPC) bilayers in the liquid crystalline phase with a number of small 
molecules including trehalose, a disaccharide of glucose, alcohols, and dimethylsulfoxide (DMSO). 
The sugar molecules tend to stabilize the structure of the bilayer as they bridge adjacent lipid head- 
groups. They do not strongly change the structure of the bilayer. Alcohols and DMSO destabilize 
the bilayer as they increase its area per molecule in the bilayer plane and decrease the order param- 
eter. Alcohols have a stronger detrimental effect than DMSO. The observables which we compare 
are the area per molecule in the plane of the bilayer, the membrane thickness, and the NMR order 
parameter of DPPC hydrocarbon tails. The area per molecule and the order parameter are very 
well correlated whereas the bilayer thickness is not necessarily correlated with them. 
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Phospholipid bilayers have been the focus of research for a long time due to their natural occurrence in cellular 
and intracellular membranes. Detailed computer simulations of phospholipid monolayers and bilayers have achieved 
a high degree of sophistication over the last years X^^^^^p-.] , 8, 9, lOj. The structures of lipid bilayers in water 



has increased significantly. More recently computational studies of mixtures of various phospholipids iiaiiaiiMlia 
and of phospholipids with cholesterol [la.[r^ 18. 19] have also been reported. However, the question how such lipid 
bilayer membranes interact with small molecules has to date not gained very much attention except for a few initial 
studies on alcohols [l^,!!!!!!, sugars [Hill, and dimethylsulfoxide (DMSO) IMil. The influence of lipophilic 
polymers on lipid bilayers has also been studied p^ . 

If we want to understand the interactions of cells with their environment, we first have to understand the interaction 
of cell membranes - or phospholipid bilayers as model systems - with molecules typical in cell environments. Cell 
membranes are the first part of the cell to come into contact with any nutrient, pathogen, or other molecule in 
the environment. So the understanding of membrane interaction with small molecules is of tremendous biological 
importance. 

In this contribution we compare the structural influences of different small molecules on model lipid membranes. 
We are focusing on molecules with high biological relevance. These include trehalose, alcohols and dimethylsulfoxide. 
Su gar molecules are obviously nutrients to living organisms and some sugars are also known to be cryo-protectants 
I2^l30j . In particular, trehalose, which is a disaccharide of glucose has been found to be very effective in this respect. 
It has been pointed out that trehalose is able to form only one internal hydrogen bond in contrast to, e.g., sucrose 
which forms two [sj- This leaves trehalose more susceptible to hydrogen bonding with lipids. Recently, it has been 
shown that the molecular mechanism underlying this cryo-protective effect is the hydrogen bonding pattern of the 



I. INTRODUCTION 
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FIG. 1: The various molecules used in our simulation and a typical bilayer configuration. 



trehalose molecules to the bilayer headgroups ||23| . The sugar can "substitute" some of the hydrogen bonds normally 
provided by water and by that stabilize the fragile bilayer arrangement. Stabilization or destabilization here mean 
that the bilayer is able to withstand harsher or less harsh conditions in presence of small molecules compared to the 
pure bilayer in water. Trehalose is experimentally known to prevent the lipid from undergoing a phase transition 
under cooling, i.e., it shifts the main phase transition temperature sig nificantly |33 |. 

DMSO is known to have some cryo-protective properties as well oa\, in addition to its anti-inflammatory, analgesic, 
muscle relaxant, mast cell stimulation and collagen dissolution properties 34] . It has been reported that DMSO 
strongly changes the overall bilayer structure by penetrating into the hydrocarbon layer to a much higher degree than 
sugar or alcohol molecules . 

It has been observed experimentally that alcohols have a destabilizing effect on model membranes [3^ |3^. It 
has also been observed that upon addition of alcohols, the lipid bilayer becomes thinner and the area per molecule 
increases. A related industrial question is the problem known as stuck fermentation in the wine industry Issf. 
A stuck fermentation means that the yeast cells do not transform all the available sugar into alcohol but stop at an 
incomplete stage. It has been proposed that the underlying mechanism is an alcohol-triggered structural transition in 
the membrane which forces trans-membrane proteins to change conformation rendering them dysfunctional [3^ . To 
date, there is no predictive means to detect stuck fermentations. The first approach to understand the mechanisms 
leading to stuck fermentation is therefore to study the influence of alcohols on lipid bilayers as the direct effect on the 
proteins is expected to be only secondary. 



II. SIMULATION DETAILS 



In this contribution we focus on a comparison between small molecules and their effect on the bilayer structure. 
The details of the simulation models are presented elsewhere '21, 23, 25] . Let us summarize the main characteristics. 
All simulations contain fully hydrated dipalmtoylphosphatidylcholine (DPPC) bilayers with 128 molecules, i.e., 64 per 
leaflet (DPPC is one of the most abundant phospholipids in animal cell membranes). Our simulations are in atomistic 
detail except for the hydrogen atoms bonded in methyl(ene) groups. These groups are collapsed into a united atom 
description centered on the respective carbon. There are at least 3655 water molecules in the system (8958 water 
molecules for some systems with alcohol). All simulations were performed at 325 K, a temperature at which DPPC 
is in the biologically most relevant fluid phase. Simulations of the pure bilayer, as well as bilayers with up to 3 wt% 
methanol or ethanol, (lipid free basis, as we do not take the lipids into account for concentration calculations) have 
been performed. Note that a weight concentration of 1.7 wt% methanol is the same molar concentration as 2.5 wt% 
ethanol. In the case of trehalose the concentrations considered were up to 3 wt%, simulations at higher concentrations 
and different temperatures have been reported elsewhere The DMSO systems contained 18.6 wt% DMSO. 

Simulations were performed under constant temperature and constant pressure conditions using the Berendsen weak- 
coupling scheme [33 . Simulations used a time-step of 1-2 fs. Figure^shows all the used molecules as well as a typical 
bilayer conflguration. 

It is worth noting that the exact simulation conditions can have a signiflcant effect on the structural properties. 
For example, changing the cutoff for the nonbonded interactions can change the area per molecule by as much as 
10 % 8]. Thus all direct comparisons have been performed under identical conditions with a cutoff of the Lennard- 
Jones interaction of 0.9 nm or 1.0 nm, respectively. The simulation models of the lipid alkyl tails used either the 
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NERD [iollZllli^ or the GROMOS force- field • n o sig nificant differences were observed between tfiese two models 
based on simulations of the pure bilayer system |3, 123 ■ This is not surprising as the two force-fields differ only 
slightly in their treatment of the hydrocarbon tails. The parameters for the head-groups are identical. However, if 
we compare our areas per molecule to the data of Pereira et al. '24] who used a different model, we see a difference 
of about 10% for the pure DPPC bilayer and for the systems with trehalose. The models for the alcohols use either 
the GROMOS forcefield ^3 or the model suggested by Miiller-Plathe 44], small differences in the density profiles 
were observed. The other observables discussed here were in all compared cases the same, a detailed comparison will 
be published elsewhere Trehalose was modeled using the OPLS model 46], DMSO using the revised model by 
Bordat et al. and water using the SPC/E model [i^ . 

Electrostatic interactions have been considered using the particle mesh Ewald technique ,49] . A few simulations have 
been performed using a reaction field technique |50l | . We compared the results from the simulations using the reaction 
field technique to the PME data for a representative subgroup and found no significant differences. Simulations did 
not contain any salt or other electrolytes. 

We focus on three key observables to characterize the structure of the bilayer systems in contact with small molecules. 
These are the area per molecule, the membrane thickness and the alkyl tail order parameter. The area per molecule 
is simply defined as the surface area in the xy plane (the membrane normal is the z direction) divided by the number 
of molecules per leafiet. The membrane thickness is defined as the distance between the peaks of highest density in a 
mass density profile, which roughly corresponds to the distance between the planes of the phosphorus atoms. 

The order parameter is defined in reference to NMR experiments by 

2 1 

— = -^Sxx + -^Syy, (1) 

Sa/3 = (3COS0Q cosO^a - 5q/3), a,f3 = x,y,z (2) 

COsOq = iaCz, (3) 

where Bz is a unit vector in the laboratory z-direction and Ba is a unit vector in the local coordinate system of the tails, 
which involves three connected carbons Ci_i,Ci, and C^+i and e = r^+i — This order parameter characterizes the 
alignment of the hydrocarbon tails with the bilayer normal. It can be used as one of the parameters (in addition to, 
e.g., the area per molecule and the thickness) to characterize the thermodynamic phase of the system. All simulations 
in this contribution are in the biologically most relevant fiuid phase. At a lower temperature the order of the tails 
increases and this triggers a phase transition to the gel phase. 

All these observables are experimentally accessible; the order parameter and the area per molecule can be measured 
by NMR. The thickness can be measured experimentally by using electron density profiles from X-ray scattering sij . 
However, especially under the influence of small molecules it is often inferred indirectly from the area per molecule 
under the assumption of volume conservation in the bilayer |35j . 



III. BILAYER STRUCTURE AND INFLUENCE OF SMALL MOLECULES 

At the temperature of 325 K we found the area per molecule to be 0.65 nm'^ for a pure DPPC bilayer which lies 
well within experimentally reported values of 0.55 nm^ to 0.71 nm^ 51, 52, 53,^]. The alkyl tail order parameters 
for the pure DPPC system, shown in Fig. |21 are also consistent with experimental data |55| and previous simulation 
studies [3- The area per molecule increases and the order parameter decreases with increasing temperature in the 
pure system Moreover, the order parameter decreases continuously from the head group to the end of the tails in 
the bilayer center. This means that the tails in the center of the bilayer are only very weakly aligned with the bilayer 
normal, especially at high alcohol and DMSO concentrations. It is conceivable that at higher alcohol concentrations 
the order may be completely lost. _At very high alcohol concentrations interdigitation between opposing monolayers 
has been experimentally observed 56]. Here, no interdigitation was observed as we study much lower concentrations 
(cf. Fig.O)). The thickness of pure bilayers shrinks as the phospholipids become less ordered. This means that in the 
z-direction - defined to be the normal to the bilayer surface - one finds a negative heat expansion coefficient. The 
overall heat expansion coefficient is positive as the increase in area overcompensates the loss in thickness. Figure 
shows the density profile of a pure DPPC bilayer. We clearly see the self-assembly of the system into a water region 
with a density very close to bulk water, an interface with a high local density as the heavy phosphorus atoms are 
located in this area, and an aliphatic region well below water density in the center. The plane of lowest density 
in the middle of the bilayer where the lipid tails from opposing leaflets meet is a symmetry plane, which we define 
as the z = plane. The interface region contains some water and the lipid head groups, which include the choline 
and phosphate groups with the density increasing up to around 1400 kg/m^. A small addition of 1.0 wt% methanol 
leads to an overall decrease of the bilayer density which suggests an increase in the area per molecule. This is indeed 
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FIG. 2: Order parameter for DPPC at 325 K in dependence of concentration of added small molecules. We average over 
both tails. Note that the data for trehalose is different from the one presented in |23l as we are presenting data at a lower 
temperature. Carbon number 1 is closest to the headgroup and carbon number 15 is the end of the tails. 



Molecule weif 


;ht % area per molec. 


[nm^] 


pure water 





0.65 


MeOH 


1 


0.69 




1.7 


0.69 




3 


0.72 


EtOH 


2.5 


0.70 


Trehalose 


3 


0.64 


DMSO 


18 


0.72 



TABLE 1: Dependence of the area per molecule at 325 K under addition of different molecules. Note that the concentrations 
are weight % lipid free basis. The Error estimates are iO.Olnm^ unless otherwise noted. 



observed (of. Table. P), as well as an increase in the bilayer thickness, leading to an overall decrease of bilayer density. 
This result actually challenges the typical experimental assumption of constant bilayer density often used to calculate 
changes in layer thickness 35, 36]. To characterize the structure of the bilayer at least two independent measurements 
out of the set thickness, area, and density are necessary. 

Experimentally it is found that low molecular weight alcohols increase the area per molecule and decrease the order. 
The experimental increase for an SOPC bilayer is of the order of AA/A « 0.1 at 5 wt% methanol. Consequently, 
the thickness of the membrane is inferred to be thinned out. Measurements using DPPC or other fully saturated 
phosphatidylcholines have not yet been reported. The only simulation study of alcohols on bilayers we are aware of |23| 
did not allow for changes in the box-size in the plane of the membrane as that study focused on local structural effects 
and dynamics, and hence it is not possible to directly compare our results to those of Feller et al. j23|. Here we find a 
linear dependence of the area expansion on alcohol concentration as well as an increase of the effect from methanol to 
ethanol which is in line with the experimentally observed increase of the area expansion with alcohol chain length [s^ . 
A more detailed analysis of the alcohol-bilayer interaction is found in Ref . . 

Bemporad et al. measured free energies of transition through a lipid bilayer for various small molecules |22l |. 
Although their numbers for the permeability are an order of magnitude too large, it is clear that methanol faces a 
barrier which is almost as high as the one for water. It is highly unlikely for the methanols to cross the bilayer. This 
is very much in agreement with our density profiles as we do not find any methanol molecules inside the layer. For 
ethanol, however, a few molecules were found inside the layer in agreement with the old Traube rule |57j | which says 
that interfacial tension drops by a constant factor as one adds methylene groups. 

The simulations with trehalose revealed only small changes to the lipid bilayer. We found a slight increase in 
the order parameter with increasing sugar concentration; this effect was earlier found to be more pronounced for 
trehalose than for its isomer sucrose |23| . We find that trehalose stabilizes and preserves the membrane under cooling 
without affecting its structure. This is understandable from the point of view of cryo-preservation and also due to 
the fact taht sugars are an abundant nutrient of cells and membrane integrity must not be strongly affected by their 
presence (see ref.|23| for additional results and discussions). As seen from the density profiles in Fig. O;, trehalose 
does not penetrate into the bilayer and accumulates at the headgroups. The sugar molecules hydrogen bond to the 
lipid headgroups, one trehalose molecule can bind to up to three lipids. 

Pereira et al. performed simulations of a DPPC bilayer in water with much higher concentrations of trehalose and 
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FIG. 3: Density profiles a) of a pure DPPC system and its constituents b) pure system in comparison witli a system containing 
1 wt % metfianol (using tfie model from Ref. |43 . We only show the density profile of the lipids here, c) system containing 
3.4 wt % trehalose, these curves are centered around the middle of the water layer {z = 3.47 nm) for clarity. Note that the 
asymmetry in the trehalose distribution is a statistical effect and not a sign of insufficient equilibration. 



found very similar effects |2J]. The area per molecule is almost unaffected even with 128 molecules of trehalose for 
128 molecules of DPPC. The same is true for the order parameter. Their simulation model is different so that the 
area per molecule cannot be compared directly, they get 0.58 nm^ for the pure DPPC system at 325 K. But it is clear 
that the effect of trehalose is very similar in two independent studies with distinctly different models. 

For the bilayer system containing DMSO, the area per molecule is increased, but not as strongly as in the presence 
of methanol or ethanol, and consequently we observe a decrease in the order parameter. A concentration of about 
18 wt% DMSO is needed to get the same increase in area per molecule as with 3 wt% methanol. Even if the molar 
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concentrations are considered for comparison, we obtain the same effect with 1.6 mol% methanol versus 5 mol% 
DMSO. 

As described earher, we see an increase in the thickness of the bilayer upon addition of methanoL We measured 
the thickness of the bilayer as the distance between the planes of highest density of the overall system. In the pure 
system (DPPC and water only) we obtained a layer thickness of 3.68 nm. Addition of 1 wt% methanol gave 4.00 nm, 
whereas 18 wt% DMSO lead to a decrease to 3.24 nm. This suggests that the change in the bilayer thickness can not 
be directly inferred from the change in the area per molecule. Again we see that trehalose has the weakest influence 
with a thickness of 3.50 nm. The errors on the thickness estimation are at most around 0.2 nm. 



IV. CONCLUSIONS 

We have found, in agreement with experimental data, that trehalose has a stabilizing effect on lipid bilayers. 
The bilayer structure is essentially unchanged from the structure of the pure system even if the temperature is 
lowered !23, 30]. However, the presence of alcohol or DMSO is contrastingly different: small concentrations of alcohol 
have a detrimental effect on the bilayer structure, the influence of DMSO being somewhat less pronounced. The area 
per molecule and bilayer thickness increase in the presence of alcohols, leading to a decrease in the overall density of 
the bilayer. The opposite effect is observed with DMSO, as the bilayer density is relatively unchanged, and the area 
per molecule increase leads to a thinning of the bilayer. 

These results suggest that the area per molecule and the alkyl tail order parameter are strongly anti-correlated, 
that is, an increase in the order parameter is accompanied by a decrease in area per molecule and vice-versa. Similar 
results have been found in a recent study on chesterol phospholipid interaction 19] . This is understood as the area 
per molecule describes the two-dimensional packing density of the headgroups and the tails. The closer the alkyl 
tails are packed, the more they are ordered. It is generally accepted that the phase transition to the gel phase - 
with a strong increase in order parameter - at lower temperatures is due to tail packing. At this moment, it is not 
clear if there is a correlation of the bilayer thickness to the degree different molecules penetrate the bilayer. However, 
even this cannot explain our observations as DMSO penetrates the bilayer better than the alcohols poll2lLl25l ]. This 
question remains open for the moment and further investigations are needed. 

In the light of our findings, that there is no direct correlation between bilayer thickness and area per molecule, we 
would recommend to use scattering experiments to obtain the layer thickness and not to rely on the assumption of 
constant bilayer volume. 

We conclude that atomistic simulations are a powerful means to study the interaction of small molecules with model 
lipid bilayers and can give us detailed insight on the local mechanisms of interactions. In order to understand the 
larger scale effects, simulations with simpler models have to be applied [sM f59ll60ll6l| . 
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List of Symbols and Abbreviations 



A 


area per molecule in the plane of the bilayer 


Sa/3 


Kronegger Delta, 5 = 1 if a = (3 otherwise 


AA 


change in the area per molecule in the plane of the bilayer 


DMSO 


dimethylsulfoxide 


DPPC 


dipalmitoyl-phosphatidylcholine 


GROMOS 


Groningen Molecular Simulation 


NERD 


Nath, Escobcdo, de Pablo, revised force-field 


NMR 


nuclear magnetic resonance 


PME 


Particle Mesh Ewald 


ScD 


NMR alkyl tail order parameter 




components of the order parameter tensor 


SOPC 


stearoyl-oleoyl-phosphatidylcholine 


z 


coordinate along the bilayer normal 


e„ 


angle between unit vector a and the bilayer normal 
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